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Abstract 

We analyze the proton decay via dimension six operators in supersymmetric SU (5)-Grand Uni- 
' fied models based on intersecting D6-brane constructions in Type IIA string theory orientifolds. 

' We include in addition to 10 1010 10 interactions also the operators arising from 5 510 10 inter- 

^ : 

, actions. We provide a detailed construction of vertex operators for any massless string excitation 



arising for arbitrary intersecting D-brane configurations in Type IIA toroidal orientifolds. In par- 
ticular, we provide explicit string vertex operators for the 10 and 5 chiral superfields and calculate 
explicitly the string theory correlation functions for above operators. In the analysis we chose the 



^ \ most symmetric configurations in order to maximize proton decay rates for the above dimension 

d ' six operators and we obtain a small enhancement relative to the field theory result. After relating 

the string proton decay rate to field theory computations the string contribution to the proton 
lifetime is t^'^ = (0.5 — 2.1) x 10^^ years, which could be up to a factor of three shorter than that 
predicted in field theory. 
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I. INTRODUCTION 

Grand unified theories (GUT's) [1] not only give a neat and aesthetic description of our 
four dimensional world but also lead to an explanation of electric charge quantization and - 
with the aid of supersymmetry - predict the value of sin^ Ow in very good agreement with the 
experimental one. Moreover GUT's lead to Baryon number violating processes; in particular 
they predict proton decay P (for a recent review on proton decay see P]). 
In supersymmetric GUT field theories 0,01 the proton decay can occur either by an exchange 
of a super heavy SUSY particle which corresponds to a decay via the dimension 5 operator 
J d?9Q^ L or by a super heavy gauge boson exchange^. The latter corresponds to a decay 

^ We forbid proton decay due to dimension four operators by introducing i?-synimetry. 
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via the dimension 6 operator J d^9 Q* L* . In the simplest supersymmetric GUT mod- 
els, proton decay mediated via dimension 5 operators dominates and recent computations 
p.adW alifetin^e fo. the p.oton, wMch i= below the present e„„tal bound. BflS, 
but [21, 110| • The fact that proton decay has not yet been observed, suggests the existence 



of some mechanism that suppresses or even forbids these dimension 5 operators 
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III, so 



121 is the most dominant one. 



that after all the proton decay via dimension six operators 
In this paper we investigate proton decay via dimension six operators in supersymmetric 
GUT models based on intersecting D6-brane constructions on type IIA string theory ori- 
entifolds. More precisely, we compute the string effects on the proton's decay into a pion 
and a positron {p 7r°e^) for supersymmetric 5'f/(5)-GUT-like models arising from inter- 
secting D6-brane constructions. In S'f/(5)-GUT's there are two different amplitudes that 
contribute to this proton decay rate: (10*1010 10) and (5*510*10), where 5 and 10 de- 
note the multiplets of the gauge group SUib). For intersecting D6-brane constructions with 

nnn n 

supersymmetric S'f/(5)-GUT's (ij, llJ, ll^^; the latter amplitude was computed in |23], by 
explicitly calculating the string amplitude contribution to 10 1010 10 operators. However, 
even after pushing all the parameters to the limit, in order to maximize the proton decay 
rate, the string contribution to it is at most comparable to the field theory one. In this 
work, we explicitly evaluate the amplitude (5*510*10) in the same class of models. 
A.in0,i;.ead of petroling the calclation in a specific model, we rather use generic 
universal features of intersecting D-brane model constructions which are relevant for deter- 
mining the proton decay rate. In general, the amplitude is sensitive to the local structure of 
the intersection and the way the D6-branes are wrapped around the compact space. Assum- 
ing that the size of the compactified volume is bigger than the string size, the latter effects 
can be neglected and the computation can be performed for a local D6-brane configuration 
where we do not need to worry about the embedding in the compact space. This approach 
allows us to make predictions about the proton decay rate in a general class of intersecting 
D6-brane orientifold models. In generic models the matter fields 10 and 5 are not located at 



^ For a review on intersecting D-brane constructions see, e.g., |16| |: for the original work on non- 
supersymmetric intersecting D-branes, see 
and also 0. For flipped SU(b) constructions see 
within Type II rational conformal field theories see: 
on Calabi Yau manifolds see 



, and chiral supersymmetric ones , see^^ \\^ 
For supersymmetric SU{5) GUT constructions 
241 and references therein. For a related study 
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the same intersection, which leads to an overall suppression of the amplitude (5*510*10). In 
this work, in order to maximize the effect, we assume the most symmetric case that all the 
matter arises at intersections that are on top of each other. Therefore, we rather compute 
an upper bound for the string contribution to the proton decay rate in these models than 
determining the complete amplitude which is model dependent. 

This paper is organized as follows. In section 2 we describe the local setup in which we work 
and derive the conditions on the intersecting angles, in order to obtain the matter fields in 
the representation 10 and 5 at the intersection, simultaneously. In section 3 we apply the 
prescription, given in appendix A, to construct the vertex operators for the matter fields. 
Section 4 is dedicated to the computation of the string scattering amplitudes, including their 
normalization. Section 5 states the results of the numerical analysis, while the details can be 
found in appendix B. In section 6 we relate the string theory results to the four- dimensional 
field theory and determine the implication of the string scattering amplitude to the proton 
lifetime. Finally in section 7 we present our conclusions. In appendix A we give a detailed 
description, how to construct properly vertex operators for strings stretched between two 
intersecting D-branes. 



II. SETUP 



We want to analyze proton decay which occurs due to dimension 6 operators in a local 
intersecting D6-brane configuration. Therefore, we have to consider scattering amplitudes of 
the form (5*510 10) and (10*1010 10), where 5 and 10 denote the multiplets of the gauge 
group SU{5). While the latter amplitude was already examined in we will determine 
the additional contribution to the proton decay arising from the amplitude (5*510*10). 
Since we shall only consider scattering arising at the local intersection, the first step is to 
derive conditions on the angles so that we have at the local intersection matter fields in 
the 5 and 10 representation, simultaneously. We will show that this condition is satisfied 
only for particular regions. For the explicit analysis we shall employ the toroidal orien- 
tifold construction and take the size of the tori larger than the inverse string tension, thus 
suppressing effects due to the world-sheet instantons. In this limit we shall calculate the 
four-point string amplitudes for the chiral superfields at the D6-brane intersections at the 
origin of the toroidal orientifold. In that sense the analysis can be applied as the leading 
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order calculation of string amplitudes for the states at the same D6-brane intersection within 
any orientifold construction. 

Let us briefly review the main properties of intersecting D6-models. Generically, one has 
a number of stacks of D6-branes {Ni denotes the number of D-branes for the i-th stack), 
which fill the four-dimensional Minkowski space and intersect each other in the internal 
space. Open string excitations located at the intersections correspond to four- dimensional 
chiral fermions transforming in the bifundamental representation {Ni, Nj), while open strings 
starting and ending at the same stack of D6-branes transform as seven- dimensional U{Ni) 
gauge bosons. In order to make contact with the real world, one has to compactify the six- 
dimensional internal space which leads to additional consistency conditions on the model 
called the RR tadpole conditions. D-branes act as sources for the Ramond-Ramond (RR - 
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charges which need to be canceled due to Gauss' law in the internal compact space 
Typically one introduces Orientifold six (06-) planes, not only because they carry nega- 
tive RR-charge, but also because they can maintain supersymmetry in the four- dimensional 
world, while the introduction of anti- D-branes would break all the supersymmetry. The ori- 
entifold action leads to image D6'-branes and open strings stretched between a D6-brane and 
its image transform as symmetric or anti-symmetric representation of U{Ni). As mentioned 
in the introduction, we rather investigate the proton decay amplitude in a local D6-brane 
configuration than in a specific model. In the following we discuss all the necessary ingredi- 
ents for this configuration to obtain a supersymmetric SU (5)-GUT like model (for the 
non-supersymmetric case see [3,13]). 

As explained above the analysis of the D-brane configuration we compactify the internal 
dimensions are on a factorizable six-torus T^. Later we assume that the compactification 
volume is larger than the string scale so that local effects dominate the amplitude and global 
ones can be neglected. This assumption also allows us to embed the local D-brane configu- 
ration, described below, into an arbitrary compactification manifold. 
The complex coordinates of the factorizable six-torus = x x are given by 

zi = + ix^ Z2 = x^ + ix'^ zi = x^ + ix^. 

In order to construct an SU (5) GUT model we shall consider very symmetric configurations 
of D6-branes. We take a stack 6 of M D6-branes oriented in the 0123468 directions that 
coincides with a stack a of 5 D6 branes along the 0123 directions and forms (supersymmetric) 



intersecting angles with stack b in the internal toroidal directions. The dimensions 0123 have 
an interpretation as a 3 + 1 dimensional intersecting brane world. Both types of D-branes 
are wrapped on the {n^ , m^) cycle of the torus. Obviously, the wrapping numbers of the 
stack b are given by 

b: {nlml){nlml){nl,ml) = (1, 0)(1, 0)(1, 0), (2.1) 

while the one from stack a can take the general form 

a: {nl,ml){nl,ml){nl,ml). (2.2) 

Given the wrapping numbers, one can compute the intersection angles which are in general 
given by (-R(,-R| denote the radii of the P'^ torus)^ 

and in our case take the simple form (since 9b = 0) 

«i = arota„(<|). (2.3) 

In order to cancel the RR-tadpoles, we must introduce 06-planes and in particular the 
orientifold action QR, where Q is the world-sheet parity and it! acts by 

This orientifold action forces us to include stacks of image D-branes. Since we chose stack 
b to lie on top of the orientifold 06-plane, it is invariant under the orientifold action: for 
M coincident branes on top of the 06-plane the QR projection leads to the gauge group 
Sp(2M). For the stack a we have to introduce an image stack a' of 5 D6- branes whose 
wrapping numbers are given by 

a' : {nl -ml) (n^, -ml) {nl -ml) . (2.4) 

Fermions that arise from strings stretched between a and a' transform in the antisymmetric 
representation of SU{5) x SU{5), due to the fact that the D-branes intersect at the origin of 
the torus. Depending on the sign of the intersection number these fermions transform as lO's 

^ Note that with this definition clockwise angles are positive and counter-clockwise negative. 
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or lO's . Fermions in the ab and ab' sector transform in the bifundamental representation 
(5, M) or (5, M)"^ again depending on the sign of the intersection number. In general, the 
intersection number for two intersecting D-branes a and b is given by 

3 

lab = 11 {nimi-mini) . (2.5) 
1=1 

Now we have all the ingredients to determine the conditions the intersection angles 6i have 
to satisfy in order to observe matter fields transforming as 5 and 10 at the intersection, 
simultaneously. Using ()2.H) . ()2.2|) . ()2.4|) and ()2.5p we obtain for the intersection numbers 

lab and laa' 

3 3 

lab = (-1)' n = (-2)' n ^'^i ■ (2.6) 

7=1 7=1 

Obviously, the sign of the intersection number depends on the sign of the wrapping numbers. 
For every angle dj (from now on we denote 6'^^ by 6j where 6'^^ is given by (j2.3|) ) we have to 
distinguish between four different cases 

• ra^, > which corresponds to an angle with < 6j < — 

IT 

• > 0, < which corresponds to an angle with — — < 9j < 

TV 

• < > which corresponds to an angle with — < 6i < 7c (2-7) 

• n^, < which corresponds to an angle with < di < — tt . 

Since we want to analyze proton decay in a supersymmetric GUT model the choice of the 
intersection angles Qi is not arbitrary; the sum has to satisfy jsol 

^1 + ^2 + ^3 = mod 27r . (2.8) 

This requirement restricts the choice of the angles. First we consider the case that the angles 
add up to and later on we also analyze the configuration where the sums of the angles are 
27r or — 27r. If the sum is equal to then one or two of the angles have to be negative. If 
only one angle is negative, let us assume without loss of generality that < 0. Since for all 
angles \Qi\ < vr, we distinguish between four different cases for which we obtain, by applying 
(12. 6p and (j2.7|) . the intersection numbers and in particular their signs 



M denotes the representation of the gauge group Sp{2M). 
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• 7„b > and > for < < f < ^2 < f - f < ^3 < 

• 46 > and < for < ^1 < f < ^2 < f - tt < ^3 < -f 

• Iab> and > for | < 6'i < tt < 62 < ^ - n < 9^ < 

• Iab> and > for < 6I1 < | ^ < $2 < tt - tt < 6^3 < -| . 

For all combinations of ^/'s that fulfill the above stated properties 9i — and one angle is 
negative) we see that the strings stretched between D-branes a and b transform as 5 instead 
of the desired 5. Therefore we do not observe a 4-point interaction of the form 5*510*10 at 
the intersection. 

Analyzing the case of two negative angles (without loss of generality we assume that 61 and 
02 are negative) we again distinguish between four different cases 

• lab<0 and < for -f < ^1 < - f < ^2 < < ^1 < f 

• lab<0 and > for -f < ^1 < - f < ^2 < f < ^3 < tt 

• lab<0 and < for -tt < ^1 < -f - f < ^2 < ^ < 63 < n 

• lab <0 and < for < 61 < - tt < 62 < ~ ^ < 63 < tt . 

Only in the region — | < 61^2 < , | < ^3 < tt we observe matter fields transforming as 5 
and 10, where strings stretched between the D-branes a and b transform as 5 and strings 
stretched between a and a' transform as 10. 

Let us now turn to the case in which the intersection angles 9i add up to 27r. Then all 
the angles are positive and we have to distinguish between three different configurations 
(without loss of generality let us assume that 9i is always bigger than |) 

• /ab < and < for | < 6*1 < tt | < 6*2 < tt < 6*3 < f 

• 4b < and 4„/ < for f < ^1 < tt < ^2 < f ^ < 63 < tt 

• lab<0 and > for ^ < 9i < tt | < 6I2 < tt f < 6I3 < tt . 

Again only in one region, | < ^1,2,3 < tt, we observe matter fields transforming as 5 and 10, 
where strings stretched between the D-branes a and b transform as 5 and strings stretched 
between a and a' transform as 10. 
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Finally, we examine the case in which the angles add up to — 27r. Here all three angles have 
to be negative and again one has to distinguish between three different cases (without loss 
of generality we assume that 9i is smaller than — |) 

• Iab> and I aa' > ior -tt < < | - tt < 6*2 < f - f < 6*3 < 

• Iab> and > for -n < 61 < ^ - | < ^2 < - tt < ^3 < -| 

• Iab> and laa' < for -vr < 6^1 < -| - tt < 6*2 < -f - tt < 6*3 < -| . 

As in the first case, the analysis shows that strings stretched between D-branes a and b 
transform as 5 under the U{5) gauge group. Therefore, at the intersection we do not have 
any matter fields transforming as 5. 

Summarizing, we determined that only for the two regions — f < 9i,2 < , ^ < 63 < % 
and ^ < ^1,2,3 < vr we have matter fields transforming as 5 and 10 at the intersection 
simultaneously. In addition to the amplitude (10*1010*10), we have for these two regions 
only, a non-suppressed contribution from (5*510 10) to the proton decay rate. In order to 
compute these two amplitudes we need the corresponding vertex operators to the states 5, 
and 10 in the respective configurations, which we determine in the next section. 

III. VERTEX OPERATORS 

For different D-brane configurations we have different vacua and therefore different vertex 
operators. Knowing the D-brane configuration we can use the prescription given in appendix 
A to obtain the vertex operator for the massless fermion in the R-sector. In this way we 
can easily determine the vertex operators for 5, arising from strings stretched between the 
stacks a and b. The vertex operator for 10 requires more effort. The simple approach just to 
replace the 9j in the 5 vertex operator by the double, 29 j only works for \9j\ < since in 
the expansion of the bosonic ()A2jl and fermionic degrees of freedom ()A3|1 the shift number 9i 
has to be in the interval [—1,1]. Therefore if 6'/ > | we need to find an expression uj which 
lies between and 1 and describes the D-brane configuration aa'. Figure ^ which shows 
the D-brane configuration for the case — ^ < 6*1,2 < 0,| < ^3 < 1. The vertex operator 

^ From now on we replace 9] by 9i/tt so that 6j e[— 1,1]. 
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FIG. 1: Intersection angles for the case — ^ < 6i <0,— ^ < 02 <0,^ <^'3< 1- 

in the (— |)-ghost picture for the massless fermion, arising from a string stretclied between 
D-branes a and b is given by (keep in mind that 61^2 are negative) 

7=1 

Now we turn to the aa' sector in which the string state transforms as 10. We see that the 
intersection angle in the third complex dimension is given by z/3 = —2 + 263. Note that the 
intersection angle is negative and lies between —1 and 0, since 63 takes a value between 
I and 1 and therefore the corresponding vertex operator for the state 10 takes the form 

3 

i;7(^) = AiOe-^(^)^"(;2) J]ai+,,(z)e<'^^+5)^^^Wei^-^W, (3.2) 

1=1 

where the angles uj are given by 

i/i = 29i V2 = 202 1/3 = -2 + 2^3 . 

Notice, that the angles uj add up to —2 so that the SUSY condition (j2.8p is satisfied. In an 
analogous way (look at figure |2I), we obtain for the other D-brane configuration 

• i<^i<l |<^2<1 |<^3<1 

For this configuration the vertex operator that creates a string stretched between a 
and b is 

J 3 
V--2(^) = A^e-^(^)5"(z) JJai_e,(^)e-'('^-^)^^(^)e^^-^(^) . (3.3) 

7=1 
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FIG. 2: Intersection angles for the case i<^i<l,i<6'2<l,^<^3<l 

The intersection angles z// are given by 

z/i = -2 + 2^1 = -2 + 2^2 z/3 = -2 + 2^3 . 

Obviously, they are all negative, so that the vertex operator which describes the mass- 
less aa'-string in the R-sector takes the form 

V^io'(^) = A^°e-^(")5"(;2)nai+,,(z)e^('^^+^)^^H^)eiA.-x(.) _ (3^4) 

1=1 

Again the angles z// add up to —2. 

In order to calculate scattering amplitudes we also need the vertex operators for 5* and 10*. 
We obtain them by replacing the spin field by the spin field with opposite chirality and at 
the same time sending the angles 9i and vi to 1 — 61 and 1 — z//, respectively (for negative 
angle we replace 9j and uj hj —1 — 61 and —1 — uj, respectively). For these two cases we 
obtain 

• -l<ei<o -|<02<o |<^3<i 

1=1 

for 5* and 

1=1 

for 10*. 
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i<01<l i<02<l i<^3<l 



for 5* and 



for 10* . 



V--}{z) = e-^(^) ~S^{z) W ae, {z) e-K-^)^^W e'^-^(^) 

7=1 



V-) (z) = A''^ e- 1 W S,iz) n (z) e-'('^^+^)^^(^) 



k-X{z) 



1=1 



(3.7) 



(3.8) 



Finally, we will discuss the Chan-Paton factors. In a setup without orientifolds strings 
transform in the bifundamental of U{N) x U{M). As already mentioned above, the intro- 
duction of orientifolds changes the transformation behavior. The full orientifold action on 
the Chan-Paton factors takes the form 



A = -InRA^InR 



where •j^r is given by 



InR 



/ Iat \ 

Itv 

1m 
V 1a/ / 

The choice of = 5 leads to the following Chan-Paton factors for the lO's 

/ Aio 0\ 



(3.9) 



A 
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Afo 



(3.10) 




\ 0/ 

where Aio is an antisymmetric 5x5 matrix. For M we choose 1 that leads to a Sp{2) gauge 
group on the D-brane b which has two components in the fundamental representation. One 
component is associated with the matter field 5 while the other corresponds to the Higgs 
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particle. Their Chan-Paton factors take the form 

/ o\ 

Ag 

\-XT / 



A 



H 



( \ 

iJ 

-E'^ 

\ / 



(3.11) 



Here Ag and H are a 5 x 1 matrices. Aio and Ag denote the usual 10- and 5-dimensional 
representations of the S'?7(5) gauge group and H is the 5 dimensional Higgs field in the 
gauge field theory. 



IV. STRING AMPLITUDE 



Having derived the vertex operators in the previous section, we have all the ingredients 
to compute the scattering amplitudes. Assuming that the compactification volume is larger 
than the string scale worldsheet instantons are suppressed and it is sufficient to compute 
just the quantum part of the amplitudes. First we will focus on (y\ y\ V^i*\^^i) and 

2 2 2 

afterwards we will compute (V^i Vj^° V^x K^i); which was already examined in 

2 2 2 2 

The amplitude V\ V^VV]\) 

2 2 2 2 

We start with the region — | < Q\ < 0,— | < 02 < 0,| < 63 < 1 and calculate the 
amplitude 

r ^ 

/ TTdz, {V\\z,)V\iz2)VTiz,)V''liz,)) , 

I 2 2 2 2 

i=l 

where the vertex operators are in the previous section. Note that all the vertex operators 
are in the (— |)-ghost pictures, which guarantees a total ghost charge of —2 on the disk. 
Plugging in the vertex operators we see that in order to calculate the amplitude we need 
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the following correlators 

4 4 4 



i<j i<j 



(4.1) 



ft ~ ~ « — - 



where Zij denotes Zi — Zj. The correlator involving the four fermionic twist fields takes an 
easy form, since we can bosonize the spin fields 
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i=l i,j=l 
i<j 

The correlator for the bosonic twist fields is more involved. Using the stress energy tensor 
methodjthe quantum part of four bosonic twist fields with two independent angles evaluates 

to 3,Q 

{ai_g{zi) crg{z2) cri_^{z3) a^{z^)) = z^^^^'^'^ z^^''^'"^ ( JLJi j r^x) , (4.3) 

with X = ^12^ and I(x) is given by 



where 



I{x) = ^ [B,{e, u) G2{x)Hi{l -x) + B2{e, u) Gi{x)H2{l 
Zn 



X 



T{l + 9-u) ' ' T{l + v-e) 

Gi(x) = 2Fi[^,l-z/,l;x] G2{x)=2F^[l-e,u,l-x] 

Hi{x) = 2Fi[e, i-v,i + e -v-x] H2{x) = 2F,[i-e,u,i-e + v- x\ 

Applying the correlators, the amplitude becomes 

A = iTr(Af A2 A^o^ kf) u^^^U2 u^^^^u^i^T,)^^^'''^ [ ^ h 

\i=\ 



X 



A [/ (-gl, 1 + Z/l, X) I [-62, 1 + U2, X) 63, 1 + Us, X) 

l\ ' (^12^34)"'^+^ (^13^24)"'*(^14^23)"'"+^ 



14 



where s, t and u are the Mandelstam variables 



The conformal Kilhng group can be used to fix three of the vertex operator positions. A 
convenient choice is 

Zi = Q Z2 = X = 1 Z4 = Zoo = OO , 

which imphes the c-ghost contribution 

(C(0)c(l)c(^oo)>=4- 

After fixing three positions, we are left with an integral over one worldsheet variable 
A = iCa Tr (Af A^ Af A^") ui-f^U2 usj^U4{27r)^S^^^ (j^ k> 



.1=1 



In order to obtain the full amplitude we need to sum over all possible orderings 

Aotai =C (rr(Af Al Af Af + Tr{Al Af Af Af)) f dx U{x) 

J —00 

C (rr(Af A^ Af Af + rr(A^ Af ^ A^^ A^^^)) j^^ dx U{x) 
C (rr(Af A^°^ A^ Af + Tr(Af A^° A^ A^)) dx U{x) , 



with 



j=i 



and 



C/(x) 



[/ (-^1, 1 + I/i, / (-^2, 1 + 1^2: X) 63, 1 + 1^3: x) 



Calculating the traces for the third term by plugging in the respective Chan-Paton factors 
immediately shows that they vanish. Explicit computation of the traces leads to the identi- 
ties Tr{Af Al Af Af) = Tr{Al Af Af Af) and Tr(Af A| Af Af) = Tr(A| Af A^o^ Af) 

15 



and thus the amphtude takes the form 

\i=i 

X (Tr(Af Af Af)K{d,, d,, 9-,) + Tr(Af A^ Af A^" ^(^1, 9„ 9^)) , (4.4) 

with 

K{9i,92,93)= [ dxU{x) T{9i,92,9^)= f dxU{x). (4.5) 

Jo J-oo 

In the field theory, the first term corresponds to proton decay via a gauge boson, while 
the second one describes the proton decay mediated via a Higgs particle, arising from the 
Yukawa interaction 10 5 5h- 
Finally we replace the z/'s by the angles 9 

ui = 29i z/2 = 292 z/3 = -2 + 2^3 

and obtain for U 

^ [I i-Oi, 1 + 2^1, x) / {-92, 1 + 292, x)I{l- 9,, -1 + 2^3, x)]'^ 

Applying the same procedure for the other sector we obtain 
• |<^i<l \<92<l \<9^<l 

The amplitude 

/ T\dz, {V\\z{)V\{z2)VX{z,)V^',{z,)) , 

/-'--'- 2 2 2 2 

i=l 

takes the form 

Atotai =2iC^ u^Yu2 usi^u, (27r)^5(^) k}j (4.7) 

X (Tr(Af A^AfAiO)ir(^^i, 02,^3) +Tr(AfA5AfAf)m,^^2,^^3)) , 
with K and T defined in ()4.5|) and given by 

3 

U{x)=x~''''-\l-x)-°'''~^Y[[I{l-9j,-l + 29i,x)]-'^ . (4.8) 
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The amplitude (l^^^r i V}^i V^,) 

2 2 2 2 

Note that in both cases, | < 61^2 < 1 , — | < 6^3 < and | < ^^1,2,3 < 1, the vertex 
operators for the matter fields transforming as 10 take the same form. Thus the compu- 
tation of the amplitude (Vj^i* Vj'^i V^^i*V}^i) is identical for both cases. We use the same 

2 222 

correlators stated above except for the one involving the bosonic twist fields, which takes a 
simpler form, since it involves only one independent angle 

\ 9(1-9) 
Zl3 Z24: \ r - i 



with 



{(ri_g{zi) ag{z2) ai_e{z3) cre{z4)) = ( ) L ^(x) (4.9) 

. 2:12 2^14 Z23 ^34 , 



^(^) = ■ } a. 2Fi[e,l-9,l,x] 2F^[9,1- 0,1,1 -x]. 
Plugging in all the correlators and fixing three vertex operator positions we obtain 

Aotai = i C'^ {Tr{Af Al' Af Af) + Tr{Af Af Af A^°)) {27.^6^'^ (fl ^ 

\i=i 

„i 3 
X Mi7'^M2 M37mM4 / dxa;-"''-^(l-x)-"'"-^ JjL-5(i + zy^ 



Finally we replace the uj by Oi and obtain 

• -\<ei<Q -\<e2<Q |<03<1 



Atotai = iC'^Tr {AfAfAfAf + Af^Af A^°^A^°^ 

X (27r)V^) \ uiYu2U3i^.u,M{e,,e2,e3) (4.10) 



.i=l 



with 



rl ™— In ™\— 1 

Mie,,92,9s)= dx— ^ , . (4.11) 

Jo L2 {1 + 29i, x)L2{l + 202, x)L-2{-l + 26-3, x) 

|<^i<l |<^2<1 |<^3<1 



A,otai = iC'^Tr[AfAl'AfAf + Af^Af A^°U^°^ 

ui'y^U2U3'y^U4M{6i,d2,63) (4.12) 
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with 



dx 



-a' s-1 



1-x) 



-a' u—1 



(2^1 - 1, x) L5 (2^2 - 1, x) L5 (2^3 - 1, x) 



(4.13) 



The (V^° V^i does not involve an Higgs exchange, since couphngs of the form 

2 2 2 2 I T 

10 10 5h are absent due to the U{1) charge conversation [lj|. 



Normalization 



In this section we determine the two undetermined constants Ca and in the string 
amphtudes computed above. We will use the fact that even in the low energy limit the 
integrals ()4.5|) . ()4.11|) and ()4.13p are convergent in the limit x — > 0, which corresponds to 
a gauge boson exchange. Factorizing the amplitude into two three point functions allows 
us to normalize it. We start with the amplitude {V^^* V^-iV}\*V}\) and turn later to 

2 2 2 2 

The amplitude {V\* V\ V}_VV}\) 

2 2 2 2 

We first examine the limit x — > and will see that even in the low energy limit the 
integral is convergent, due to the special kinematics of this problem. 

Limit X ^ 

As X — the hypergeometric functions behave like 

F(a, 6, 1, x) ^ 1 F(a, 6, a + 6, 1 - x) ^ T^r^^ ^ 

r(a)r(6) \ X J 

with 

In K(a, h) = 2^(1) - ip{a) - ip{b) . 
Applying ()4.14|) / takes the form 

lim/(^,.,x) = ^lnf^') , 
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where In 5{9^ v) is given by 



\n5{e, u) = 2^(1) - - -9)- Up{u) - 1^(1 - u) . 

Therefore even for s = t = we obtain for the integral ()4.4|) a convergent expression in the 
hmit X ^ 



3/2 



X 



\n[l/x\ 



-3/2 



(4.15) 



That allows us to normalize the amplitude by factorizing the amplitude in the limit x — 0, 
where it reduces to a product of two three-point functions 

A,{k^, k,, k„ k,) = -j ^ . (4.16) 

The unusual factor of | is introduced to take into account the doubling in the Chan-Paton 
factors. 

The three-point amplitudes describe the exchange of a gauge boson and are given by 



A^{k,,k2,k3) = igDe (27r)45W ( E ^0 ^i7^«2Tr(Af A^) . 



(4.17) 



, i=l 



A. 



(4.18) 



Here fi corresponds to the polarization and A^ denote the Chan Paton factors of the gauge 
boson. The latter takes the form 

/ Aa 0\ 
A„ 

\ 0/ 

where the A^'s are the gauge bosons of f/(5) which satisfy Tr(Aa A;,) = ^Sab- The intermediate 
state is a massless a — a string, which is a gauge boson, that can carry arbitrary momentum 
p along the directions of the D-brane a orthogonal to the intersection. In these directions 
we have to integrate over 

d^q [ dxx"'«'-"'^-i = 7r3/2 (a')"3/2 f ^^^-^'s-i [1^(1/3;)] -3/2 
Jo Jo 

which tells us that the replacement, going from effective field theory in four dimensions to 
the form of the string integrand near a; = is no longer 

-^a' [ dxx-"'^-^ , 
s Jo 
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but 



(f - s ' ' Jo 



Performing the integral on the right hand side of ()4.16|) and using the replacement ()4.19|) 
we obtain 

2 5/2 / 4 \ r 

i^^^Tr{Af AlAf Al')u,Yu2U3l,u,6^'^ (^^0 J^dxx-'^' [\n{l/x)]-'/\ (4.20) 
This needs to be the same as ()4.4p in the limit x — > 
2iCATr(Af A^A^°^Af)Mi7'^M2M37MW4(2vr)^7r=^/2^(^) (^Ylk}j dxx"^"^ [ln(l/x)]-=^/2 ^ 



which leads us with gf)^ = (27r to the normalization constant Ca 

CA = \gsa . (4.21) 

For the second amplitude one obtains, following the same procedure, the same normalization 
constant. 



The amplitude {VX V}_'{ VXV}_'{) 

2 2 2 2 

Note that the amplitude is invariant under the exchange of x and 1 — x if one si- 
multaneously interchanges s and u. Therefore we obtain similar limits for x — and x — 1. 
That is not too surprising taking into account that we expect an exchange of a gauge boson 
in both limits. 

Limit X ^ and x 1 

Using ()4.14|) and taking the low energy limit s, t ^ we get for x — > 

3/2 [ ^ in[l/x]-3/2 (4.22) 



X 



and a similar result for x ^ 1 



^ j' _ ^)]-3/2 _ (4 23) 

Following the same procedure as in the case of the amplitude {V\ V\ l/^if ) we obtain 

2 2 2 2 

for normalization constant Ca 

C'^ = 7rgsa'. (4.24) 
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V. NUMERICAL ANALYSIS 



We want to compute the contribution of the amphtude which arises from the four-Fermi 
interaction in the low energy effective theory. That means that we take the low energy limit 
and subtract the s, t and u poles, if present. It turns out that the amplitudes are divergent 
only in the limit x —oo. As derived in appendix B there is no massless exchange in the 
ti-channel. The s-channel requires more explanation, since in general we expect a massless 
gauge boson exchange, which leads to an undesired s-pole. We saw that the integral does 
not diverge at the s-pole, since we neglected global effects coming from the internal space. 
Locally, the internal dimensions look like a flat space with infinite volume which leads to a 
vanishing gauge coupling in four dimensions 

gYA^ ~ 77- , (5.1) 

here Vint denotes the internal volume and gvM is the gauge coupling in four dimensions. 
Thus even if we observe a gauge boson exchange, we do not see an s-pole in our effective 
low energy theory. In the limit x —00, which corresponds to a t-pole, the integral is 
divergent and in order to obtain the four-Fermi interaction we have to subtract this pole. A 
detailed discussion of the numerical analysis of the integrals K,T and M in the amplitudes 
(I4.4|) . (j4.7j) . (j4.1(Jj) and (j4.12|) can be found in appendix B, where for simplification we set 

Oi = 02 = 6. Table HI shows the contribution M for the string amplitude {V^i* V^i V^\*V^\) 
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-i < 01 < -i < 02 < i < 03 < 1 


i<0i<l 1<02<1 i<03<l 





K 


T 


M 





K 


T 


M 


-.40 


6.5 


5.4 


10.3 


.505 


1.5 


1.5 


2.5 


-.42 


5.7 


5.1 


9.4 


.51 


2.0 


2.1 


3.5 


-.44 


4.9 


4.6 


8.3 


.52 


2.9 


2.9 


4.9 


-.46 


4.0 


4.0 


6.9 


.54 


4.0 


4.0 


6.9 


-.48 


2.9 


2.9 


4.9 


.56 


4.9 


4.6 


8.3 


-.49 


2.0 


2.1 


3.5 


.58 


5.7 


5.1 


9.4 


-.495 


1.5 


1.5 


2.5 


.60 


6.5 


5.4 


10.3 



TABLE I: Contribution to T and M for different angles 
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and the contributions K and T arise from (V^^i* V^^i Vj^°*Vj^'l) for different angles 9. For 

2 2 2 2 

9 = —1/3 and 9 = 2/3 we observe a second massless fermion which indicates that we now 
have N = 2 supersymmetry. Since our world is chiral we choose 9 in the ranges, given in 
table □ 

Note also, that going from the first sector — |<^^i<0,— i<^^2<0,i<^^3<lto the 
second one |<6'i<l,i<^^2<l,|<^3<l and replacing 6* by 1 — 9, simultaneously 
leads to the same results for K, T and M. This is not too surprising, since the respective 
vertex operators correspond to the same states if you interchange 9 with 1 — 9. 



VI. COMPARISON TO FOUR-DIMENSIONAL FIELD THEORY 



In this section we want to compare the amplitude obtained due to massive string states 
in string theory with the amplitude on the field theory side. Therefore, we would like to 
replace all the string theory parameters such as the string coupling gg or the gauge coupling 
gD6 by appropriate expressions using quantities about which^we have some knowledge of, 
such as Mqut and acuT- We follow closely the analysis of 
The action for the gauge fields living on the Z}6-branes is 



where the Fij is the Yang-Mills field strength and Tr denotes the trace in the fundamental 
representation of U{N). After compactification on R x Q the action becomes 

where Vq is the volume of Q. Keeping in mind the usual convention Tr{QaQb) = ^^ab we 
finally obtain for the action 



d^xTiF.F'^, 



Vo 



89l 

On the other hand, the GUT action is given by 

1 



d^xTiFijF'^. (6.1) 



^9gut 



d'^xTiFijF'^ , (6.2) 



where Qgut is the GUT coupling. Comparing ()6.1|) and ()6.2p . along with gjy^ = {211)^ g. 



/3/2 



and a^jj^ = g^^^ / (Air) , leads to the identification 

2/3 

\i2n)^9. 



a'=( ^f^^ ) . (6.3) 
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The volume Vq enters into the running of the SU{3) x SU{2) x U{1) gauge couphng from high 

— 1/3 

energies to low energies. Approximately, one can say that Vq plays the role of the mass 
scale unification Mqut in four dimensions. In order to obtain the exact relation between 
them one needs to compute the one loop threshold correction to the gauge coupling, which 
was done for M-theory on a manifold of G2 holonomy js^^ 

IVlQlJrp 



where L{Q) is a topological invariant, the Ray-Singer torsion. In [22| it is argued that this 
relation holds true in Type IIA string theory and thus we finally obtain 



We would like to replace all the string parameters in the amplitudes fj4.4|) and ()4.7|) in terms 
of four dimensional field theory quantities. Unfortunately, equation ()6.5p still includes two 
string parameters L{Q) and Qs- The Ray-Singer torsion L{Q) depends crucially on the 
compact space and takes for simple lens spaces values around 8 35|- In order to neglect 
higher order loop amplitudes the string coupling Qs is better smaller than 1. On the other 
hand we are interested in the largest possible contribution to the enhancement and set 
therefore Qg approximately to 1. 



Field theory amplitude 



After relating the string parameters to four dimensional field theory constants, of 
which we have some experimental knowledge, we now recall the analysis of proton decay in 
the S'f/(5) GUT modeF. This treatment closely follows P|. 

The kinetic energy for an SU{b) gauge theory, involving the gauge field A, the fermionic 
field 'j/'g, which transforms as 5, and the fermionic field ipio transforming as 10 under the 
SU{b) takes the form 

T = j^TriF'iA)) + i^p-,rD,^-5 + iV'io7'^I^>io (6.6) 

^ An explicit computation for the one loop threshold correction in type IIA string theory was performed in 

[3^ , which leads in the limit — > 1 to an equivalent relation. 
^ As done usually we neglect because of the weakness of the Yukawa couplings to light fermions the Higgs 

mediated Proton decay. 
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with 



and 



ab ^ 9 GUT 



D' ih "-^ 



By explicitly using the antisymmetry of ijjio, the latter can be simplified to 



10 



V2 



The gauge field A can be displayed as a 5 x 5 matrix 



A.. 



^1 u M u \ 



1 /^a \a 



^2 ^, 



X. 



c 



c 



X2, 



3fi 



w- 



( 



+ 



5„ 



-2 



v 



3/ 



\ ^l/i -f2pi -f3|j 

where the A" are the Gell-Mann matrices, the denote the gluon fields of ^[/(S) and VT^, 
ly^, -B^ are the bosons of the SU{2) x The X and F are the new gauge bosons 

that are contained in SU (5) and do not occur in the standard model. The exchange of these 
new gauge bosons leads to Baryon-Lepton number violating processes and therefore allows 
proton decay. 

To make contact to the standard model the SU (5) needs to be broken, which will be achieved 
by giving the Higgs field, which transforms under the 24-dimensional adjoint representation 
of SU (5) an expectation value. This generates a mass Mx of order of 10^^ Gev for the gauge 
bosons X and Y. 

From ()6.6|) one can easily deduce the effective four-Fermi interactions which lead to proton 
decay. Ignoring mixing effects as well as second and third families one obtains for the 

2 



where the first factor arises from a 10* 10 10* 10 interaction and the second factor from a 
5* 5 10* 10 interaction . 
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Comparison 



This result fl(i.7j) we want to compare with the string theory contribution. In order 
to do that we turn on Wilson lines, that break the SU{5) gauge group into the standard 
model ones. Assuming such a mechanism of symmetry breaking exist we compute the 
traces of ()4.4|) and ()4.10|) only for entries which lead to proton decay. One obtains for ()4.4|) 
and (jHH) 

= i (2vr) V^) ng^a'^e^f.^u'l^rul) {e^n^.di) {K{e) + T{e)) (6.8) 



.1=1 

and for (jOTHl and 

Comparing the string theory proton decay rate with the one from four dimensional gauge 
theory one obtains 

Most recent calculations 3 for the proton decay mediated v. gauge bosons m an .^(5)- 
GUT model gave the hfetime t^'^ in terms of gauge boson mass Mx and a^^j, 

0.04 \^ f ' 



GUT 



^ 1.6 X loVa. ( ^ [^) . (6.11) 



This leads with the values Mx = M^^,^ = 2 x 10 °GeV and a^y-j, = 0.04 to a proton 
lifetime of 1.6 x lO^^years. The present lower bound on the proton lifetime for p —>■ e'^ 
is 1.6 X lO^^years [37] and even the next generation proton decay experiments, based on 
underground water Cherenkov detectors will reach a lower bound not larger than lO^^years 
Therefore in the near future, unless there is an enhancement to the proton decay 
amplitude, we will not observe the proton decay via gauge boson exchange. Using ()6.1()|1 
and 1)6.1111 the proton lifetime in the considered type IIA string models is 



ST '^fi 542 / 0.04 \ / M^„^ 
^ 1.6 X lO'^ears jj. 7^KT1&^ 
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where ^ ' (Q)9s — ((^^+^) +'^^ ) jg the string enhancement factor. Note that in (j6.12j) the heavy 
gauge boson mass Mx, which is model dependent, is absent and the proton hfetime depends 
only on M^^^, . We also observe an anomalous power of a^jjj, in ()6.12|) indicating the stringy 
nature of the enhancement. 

Let us examine the enhancement factor ^ ^ (Q)9b — ((^^+^) As already mentioned earlier 

the Ray-Singer torsion is around 8 for lens spaces with small fundamental group. The 
string couphng takes values between and 1, but in order to obtain the largest possible 
enhancement to the proton decay amplitude we assume it is approximately 1. Table H] shows 
that M ranges between 5 and 10, while K + T ^ 1.2 x M, leading with the numerical 
four-dimensional SU{5) supersymmetric values M^^^ = 2 x lO^^GeV and a^^^^ = 0.04 to a 
proton lifetime = (0.5 — 2.1) x lO^^years. We see that although there is in addition to 
the contribution to the four-Fermi interaction which in field theory are due to gauge boson 
exchange, there is also a contribution due to terms that in field theory arise from Higgs 
particle exchange, the total string contribution is not large enough to lead to a considerable 
enhancement in the proton decay rate. 

The dimension six operators 5*510 10 have in contrast to the operators 10*1010*10 a 
second proton decay mode; they lead in addition to the decay mode p ^ tt^ also to 
p Ti^u. Plugging in the respective entries in ()4.4|1 leading to the mode p vr+z/ one 
obtains 

41T = i(2vr)^5W (^^k}j Trg^a'^e^.^ul-^rd^) (u^^.di) {K{9)+T{9)) . (6.13) 
Within the field theory the effective interaction 

2 

= |fT(^«/^7^L"7'^4) {^Rl,d%) , (6.14) 
the ratio between the proton decay rates is given by 



For this decay mode the string enhancement to the proton decay rate is even smaller than 
for the mode p — > 7r°e+ due to the absence of the 10*1010*10 interaction term. For the 
same choice of parameter as above (in addition we assume that Mx = Mgut) the ratio 
f|6.15|) takes values between 0.2 and 0.8. 
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VII. CONCLUSIONS 



In this paper we computed the local, string contribution to the proton decay rate for 
supersymmetric SU(5) GUT's based on intersecting D6-brane constructions in Type IIA 
string theory orientifolds by explicitly calculating the string amplitude contribution to the 
dimension six operators. If the compactification volume is larger than the string scale, world- 
sheet instanton effects are negligible and the local contribution is the dominant one. In the 
computation presented, we assumed that the matter fields 5 and 10 are located at the same 
intersections on top of each other, and thus the leading string amplitude contributions have 
no suppressions from area factors. In this case the amplitudes give the largest possible con- 
tribution to the proton decay rate. In contrast to the authors j27|, who only considered the 
amplitude (10*1010*10), we also included the explicit calculation of the string amplitude 
for 5*510*10 operators. 

As a by-product we explicitly constructed the vertex operators for any massless string ex- 
citation at supersymmetric D-brane intersections arising in Type IIA toroidal orientifolds. 
Specifically, by employing explicit string vertex operators for the 10 and 5 chiral super- 
fields, we calculated explicitly string theory amplitudes contributing to the proton decay 
via dimension six operators. In the analysis we chose the most symmetric configurations in 
order to maximize proton decay rates for the above dimension six operators and we obtain 
a small enhancement relative to the field theory result. In contrast to the string amplitude 
(10 1010 10), where only the gauge boson exchange contributes to the proton decay rate 
for the amplitude (5 510 10) there is an additional contribution corresponding to the pro- 
ton decay mediated via Higgs particle. 

After relating the string theory result to the field theory computations we obtain for the 
proton lifetime in type IIA string theory models 



54^ (_04_Y' f^aur_\ 

T ((K + TY + 4M2) Wut) \2Q^'GeV) 



4 



t'/ ^ 1.6 X 10^Var. ^.,,„, . , [-^] ( ) , (7.1) 



L^/\Q)g7'{{K + TY + m^) 

which has an anomalous power of acuT indicating the string effects. The string enhance- 
ment factor depends on the Ray-Singer torsion, the string coupling gg and the numerical 
quantities M, K and T. Here the quantity M corresponds to the contribution arising 
from the string amplitude (10 1010 10), while the sum K + T originates from the string 
amplitude (5*510*10), where K is the contribution due to the gauge boson exchange 
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and T describes the contribution due to the Higgs particle exchange. Choosing common 
values for L{Q), assuming that the string coupling Qs is approximately 1 and plugging in 
the computed numerical quantities K, M and T (see table H} the proton lifetime (j7.1|) is 
Tp"^ = (0.5 — 2.1) X lO^^years, and could lead up to a factor of three shorter lifetime than 
that predicted in field theory. 



Acknowledgements 

We would like to thank Carlo Angelantonj, Andre Brown, Peng Gao, Paul Langacker, 
Tao Liu and Stephan Stieberger for useful discussions. The work is supported by an DOE 
grant DE-FG03-95ER40917 and by the Fay R. and Eugene L. Langberg Chair. 



28 



APPENDIX A: VERTEX OPERATORS FOR INTERSECTING D-BRANES 



This appendix discusses the vertex operators of bosonic and fermionic string states arising 
in intersecting D-branes based on the example of intersecting D6-branes. In the following 
we will consider D6-branes in flat, non-compact Minkowski space that fill out the first 
four dimensions (our actual spacetime) and intersect in the 3rd, 4th and 5th complex plane. 
Strings that are stretched between these D-branes have to satisfy special boundary conditions 
in the internal dimensions which leads to non-integer mode expansions for the degrees of 
freedom. In the vertex operators for the corresponding string configuration on introduces 
bosonic and fermionic twist fields to take into account these non-integer mode excitations. 
These twist fields depend crucially on the choice of intersecting angles. In this section we 
will present a instruction to construct the vertex operators arising from strings stretched 
between intersecting D-branes in the NS-sector as well as in the R-sector. 
As a first step we deduce the mode expansions for the bosonic and fermionic degrees of 
freedom. We start with the NS-sector, where strings stretched between the intersecting D- 
branes correspond to massive scalars in the four-dimensional space-time. After deriving the 
mode expansions we quantize the string, impose the condition for physical states, and obtain 
the mass formula. Later we will also deal with strings in the R-sector and show that in this 
sector we always have a massless fermion, independent of the choice of the intersection angles, 
while in the NS-sector the scalars become massless only for particular choices of angles that 
match with the supersymmetry condition. To get an idea of how the vertex operators look 
like, in particular in the internal dimensions, we examine the operator product expansions 
(OPE's) of the bosonic and fermionic fields with specific string excitations. These OPE's 
show the same behavior as the OPE's of the twist fields in orbifold theories Q]. Therefore 
the vertex operators for strings stretched between intersecting D-branes will involve bosonic 
and fermionic twist fields, ag and sq in the internal dimensions. The exact knowledge of the 
OPE's of the bosonic and fermionic fields with the string states allows us to write the vertex 
operators for the string states in arbitrary intersecting D-brane configurations. 
An open string stretched between two D-branes at an angle nOi has to fulfill the boundary 
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conditions 



4l| 



9,XP(r,0) = = X*'+^(r,0) 
(9^Xp(r, tt) + tan (nOi) (9^Xp+i(r, tt) = 
Xp4.i(r, tt) — tan (tt^^/) Xp(r, vr) = . 



(Al) 



Given these boundary conditions, we can deduce the mode expansion for (we use complex 
coordinates = X^^+^ + iX^^+^) to 



a, 



n+0i ^-"-fj 



^ for / = 1, 2, 3 . (A2) 



Z\zrz) = Y. 



n+di ^-ri-di _|_ ^ '^n-ei --n+di 



a. 



■li It, 
Upon quantization the only nonvanishing commutator is 



World-sheet supersymmetry 



leads to the same modding for the complexified worldsheet fermions (here we already used 
the doubling trick) 



(A3) 



Notice that we consider the NS-sector where the fermions are half integer modded. The only 
nonvanishing anti-commutator is given by 

For positve 6j [0 < 6j < 1) the vacuum in the internal dimensions is defined by 

1 



«LeJO) = m>l ^;^,J0) = r> 



0) = m>0 \0) = r>-. 



(A4) 
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The physical state constraint requires annihilation with all the positive modes of the Virasoro 
generators L„, in particular with Lq, which takes the form 

3 



fj.=0 \ntZ ntZ 



n 



3 



. (A5) 

1=1 KmeZ meZ ) 

Here and ip!^ denote the excitations in space-time and eo is the zero point energy. Using 
the fact that the zero mode represents the momentum of the string we manipulate 
equation ()A5|) and obtain a mass formula for the open string in the twisted sector 



11=0 \neZri neZ 
1=1 [meZ meZ 



(A6) 



The zero point energy can be computed from the (^-function regularization, as we demon- 
strate in the following (for one internal dimension only) 

-1/2 

4= [<^-m+ej,am-ej]+ ^ {r - Oi) {^p_r+ej,A-ej} 

m=—oo m=—oo (^J^7^ 

= C[-i,ej]-c[-i,i/2 + 0i] = -^ + ^ej . 

To get an expression for the vertex operators we need to determine the OPE's of and 
with some particular excitations. First we examine the vacuum state | 0) 

1 

OO 2 
r=—oo r=— OO 

where t/(0) denotes the excited twist field at the intersection. Similarly we obtain for 
Using the same procedure, the OPE of \1' and ^ with the state ip^i^g^ \ 0) is 
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Considering a negative angle 6i {—1 < 6i < 0) leads to a different definition of the vacuum 

1 



«^-eJO) = m>0 ^l^^\0) = r>- 

and the zero point energy, calculated in the same way as above, takes the form 



(A8) 



(A9) 



(keep in mind, that the angle dj is negative). Again we examine the OPE's of some special 
physical states with the fermionic fields "^{z) and ^{z). For | 0) we get 

and for ip-i-dj I 0) 

^'(^) ^-1-., 1 0) - z'^'' t,(o) ^\z) 1 0) z-'-'^' m . 

Before formulating the vertex operators for particular states we also need the OPE's with 
the bosonic fields 

dZ'{z) I 0) ^ z-'^^-''^ Ti{0) dZ\z) I 0) ^ z-'' r,(0) 

d Z\z) I 0) ^ z-'' r^O) d Z\z) I 0) ^ z-^^-''^ r^O) . 

For negative angle, we replace 6i hy ai = 1 + Oi. 

Now we can start to construct the vertex operators for the respective states. First we 
consider the state x = i^-i^e^ I 0)' where 9i,92 are negative and 63 is positive, which means 
that the string starts at D-brane a and ends at D-brane b (see figure Q)^. The mass of this 
state is given by 

3 

M' = -l-le,-^-9,+ ^-9s + l-es = -Y,9, . 

The scalar x becomes massless when the sum of the angles adds up to zero. This is in 
agreement with the supersymmetry condition. The corresponding vertex operator in the 

(-l)-ghost picture takes the form 

2 

V^\z) = e-'^(^) JJae,(;2)e'^^^^(^Vi+e3(z)e^(i+''3)^3(")e''=-^(") , (AlO) 
1=1 



Recall that we count counter-clockwise angles positive. 
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where the if/'s denote the bosonized worldsheet fermion Notice that in the case of 
supersymmetry, when the state becomes massless (/c^ = 0), the conformal weight of the 
vertex operator adds up, as required, to one. 

The corresponding complex conjugate state x* is represented by the same excitation as above 
but oriented from brane h to brane a. That means that the intersection angles 9\ = —Oj 
take the opposite sign as before and therefore the vertex operator is given by 

2 

V-\z) =e-*(") JJai_e,(2)e-'^^^^(")a_e3(;z)e-'(i+^^)^^(^)e'*=-^(") , (All) 
1=1 

Let us take a closer look at the vertex operators in the case of supersymmetry, when they 
carry a = 2 world sheet charge H = Yl'i=i Hi- The chiral superfield x has N = 2 world 
sheet charge +1, while the charge for the complex conjugate partner x* is -1. 
Next, we examine the state x* = V'-i+ei '^-i+e2 "^-l+Sal 0); where < < 1 for all /. 
Again the string is oriented from brane a to brane b (see figure |21). Why we denote the state 
by X* rather than x becomes clear later. The mass of x* is given by 

M' = -l + lt0--t(-l+0-)='-'2t0' 

7=1 7=1 ^ ^ 7=1 

and becomes massless, when the sum of the angles is equal to two, again in agreement with 
the supersymmetry condition. The vertex operator in the (-l)-ghost picture corresponding 
to this state takes the form 

V^:'\z) = e-'^(^) n ae.iz) e'(^-^)^^(^) e^'^'^^^) , (A13) 

7=1 

and as above the requirement that the vertex operator has conformal weight one is satisfied. 
The corresponding complex conjugated state x is stretched from brane b to brane a and the 
intersection angles 6'j = —6j are all negative. Therefore the vertex operator is given by 

3 
7=1 

A look at the N=2 world sheet charge in the case of supersymmetry (^/^^ Oj = 2) explains 
the notation since x* carries charge -1 while x carries +1. 

We now turn to the Ramond sector, in which the string excitations between two intersecting 
D-branes correspond to space-time fermions. The mode expansion for the fermionic degrees 
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of freedom takes the same form as for the Neveu-Schwarz (NS)- sector, but now we sum over 
integers instead of half integers 

n n 

Nothing changes for the bosonic world sheet fields Z{z,z) and Z{z,z). The vacuum is 
defined by (0 < < 1) 

«^-^jO) = m>l V'.'-ejO) = r>l 

(A16) 

a^+ejO) = m>0 V'.V|0) = ^>0- 

With this definition the zero point energy is independently of the choice of angles given by 

eo = , (A17) 

and therefore we always have a massless fermion in space time. While the mass of the vacuum 
is independent on the angles the vertex operator for the vacuum |0) depends crucially on 
the choice of angles. Let us therefore examine the OPE's of worldsheet fermions^ with |0) 
for the two different situations that we have positive and negative intersecting angles. We 
obtain for < 6'/ < 1 

^^(z) I 0) ^ z-^+^' ti{0) ¥{z) I 0) ^ z^-''' ti{0) . 



For negative angles we must change the definition of the vacuum to 

a^_,jO) = m>0 V'^^_,J0) = r>0 

«^+ejO) = m>l V'.V|0) = 



(A18) 



The zero point energy is still zero. But now we obtain different OPE's for the vacuum | 0) 

■^'{z) I 0) ^ z"2+'' ti{0) ¥{z) I 0) ^ z-'^-^' ti{0) . 

As before for the NS-sector we present for particular states the vertex operators. The first 

state we consider is the vacuum state x = I 0)) whose mass is independent of the choice of 
angles equal to zero. Assuming, that the intersecting angles 9i, 62 in the first two internal 
dimensions are positive and 6^ negative, the vertex operator takes the form 

V-'^z) = e-'^^^^ S^{z)f[ag^{z)e<''~^)'''^'^ a,^^^^^^ , (A19) 

7=1 



^ The OPE with bosonic world-sheet fields is the same as before for the NS-sector. 
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, ]_<xr I lay- 
where 5"^ = e 2 1 2 2 denotes the spin field with positive chirahty^'^. As for the NS- 

sector the corresponding vertex operator for the complex conjugated state x* is simply 

given by orientation reversal, so that the intersection angles are 9'j = —9i. Thus the vertex 

operator in (— |)-ghost picture has the form 

V;}{z) ^ e--^^'^ ~S^{z)Y[a^_eA^)e<''-^)'''^'^ a^e,^^^ , (A20) 

1=1 

.lor -J- lor 

where Sq. = e ^ ^ 2 2 represents the spin field with opposite chirahty as 5"^. Notice that 
independent of the choice of angles the vertex operator has as expected conformal weight 
one. As expected, in case of supersymmetry (^7=3 9i = 0) the vertex operators x X* 
carry N=2 world sheet charge — | and |, respectively. 

Finally let us assume that all the intersecting angles 9i are positive. In that case the vertex 

operator for the vacuum state x* takes a very symmetric form 

^ 3 
V~~'{z) = e-^(^) S,,{z) II ae.iz) eK^^-^)^^(-) e'^'^W . (A21) 

1=1 

For a similar reason as in the NS-sector we call this vacuum state rather x* than x, since in 
case of supersymmetry (X]7=i = 2) it carries | N=2 world sheet charge. Following the 
procedure described above we obtain for x 

^ 3 

VP{z) = e-2(^) JJai_e,(^)e-'('^-^)^^(")e''=-^(^) . (A22) 

1=1 

One can easily check that in case supersymmetry the vertex operator carries as expected 
N=2 world sheet charge H = 



gKi,2 g^j.g i^jjg bosonized world sheet fermions where a denotes the four dimensfonal complexified indices. 
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APPENDIX B: NUMERICAL ANALYSIS 



Before we extract the low energy limit of the amplitudes, given above, let us take a look 
at three different limits, namely x — > 0, x — > 1 and x — > — oo. The first one corresponds 
in the field theory to a gauge boson exchange, while the latter one corresponds to a Higgs 
boson exchange. In the limit x ^ 1 the type of the exchange particle depends on which 
amplitude we examine; it is either a massive particle, for {V\ V^V^VV^l) or agam a 

2 2 2 2 

gauge boson for {V}^,* V^^l VVVl^i)- We start with {V\* V\ VVVl'i) and turn later to 

2222 2222 

2 2 2 2 



X 



The limit a; ^ was already explored in section 4 in order to normalize the ampli- 
tude. Here we just state the result for the case that 6i = 62 — 



■\<9i<Q 



i<^2<0 i<^3<l 



Jo 



dx 



X 



In 



where In S{9, v) is given by 



X 



In 



5(1 + 2^,-1-4^) 



X 



In 5{d, v) = 2^(1) - - 1^(1 - ^) - iv^(i.) - 1^(1 - ly) . 



(Bl) 



(B2) 



• i<^l<l i<^2<l i<^3<l 



~ TT' 



3/2 



dx 

X 



In 



5{-e,-l + 29) 



X 



In 



5(-l + 2^,3-4^) 



X 



(B3) 



with the same S{9, u) as above. 
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Using the properties of the Hypergeometric function, in particular the transformation law 

7^/7 \ T(c)T(c — a — b) ^, , , . 

2Fi(a, 6, c; x) = — — 2Fi{a,b,a + b - c + 1;1 - x) 

r(c — a) r(c — b) 

+ (1 - zy-'^-' ^^^y^l^^ + ^ ~ "^ Fi(c - a, c - 6, c - a - b - 1; 1 - x) 

and the limit 

lim 2Fi{a,b,c; x) = 1 , 

we obtain 

I r(l-a)r(b)r(l+a-fe) xfc_a , 

T-1/ 7 \ j r(a)r(i-b)r(b~a) -^^ u^u 

r{i~a)r{b)r{a-b) ^ ^)'' * a > & 



lim-/-i(a,6,a:)^<' r(a)r(i-.)r(6~a) - - - _ ^^^^ 

2n r(a)r(i-b)r(i-a.+fe) 



In this limit we do not obtain an integer mode, which tells us that the exchange particle is 
massive. The mass depends on the choice of angles, as we will show based on our first case 
(— I < ^1 < , — I < < , 1 < < 1). Let us assume that the two angles 9i and 02 are 
equal 

6^ = 62 = 6 ^6;^ = -26 . (B5) 
In the limit x ^ 1 the amplitude ()4.4j) takes the form 

r(i + g)r(i + 2g)r(-3g) / r(i + 2g)r(2 + 4g)r(-i-6g) 
r(-0)r(-20)r(i + 3^) y r(-20)r(-i -4^)r(2 + 6^) 

for > -1/3 and 



J (1 - a;)-"'"+^+6'' (B6) 



Ti-9)Ti-29)Ti2 + 39) / r(-2^)r(-l - 40)1(3 + 69) f _ ^^_,,„_3-6e 



r(i + 9)T{i + 20)r(-i - 39) y r(i + 2^)1(2 + 40)r(-2 - q9) 

(B7) 

for 9 < —1/3. In the low energy limit ()B6|) and ()B7|) are proportional to 

a'u — a'M'^ 



37 



where M denotes the mass of the exchanged particle and is given by 

a'M^ = { \ (B9) 




which becomes massless for 6 = 1/3. For this choice of angle we observe N = 2 supersym- 
metry in the Minkowski-space. Since we focus on models with N=l chiral fermion sector, 
only, we do not take this limit. For our second amplitude ()4.7|1 we also observe a massive 
particle exchange in this limit 



r(e)r(-l + 29)T{3 - 39) /r(-l + 2e)T{-2 + 4^)r(5 - 66) 



r(i - e)T{2 - 2e)T{-2 + se) y r(2 - 2^)1(3 - Ae)T{-A + ee) 



for e > 2/3 and 



r(i - e)T{2 - 2e)T{-i + 39) / r(2 - 2^)1(3 - 4^)r(-3 + 6^) 



r(^)r(-i + 2^)r(2 - 3^) V r(-i + 2^)r(-2 + 4^)r(4 - 6^) 



(BIO) 



(Bll) 



for 6 < 2/3. In our effective low energy theory we integrate out all massive states, so that 
the part of the amplitude arising from these string massive state exchanges contribute to 
the four-Fermi contact term. 

X — > —00 

At last let us examine the limit x — oo. As mentioned earlier the second terms 
of ()4.4|) and ()4.7p give the contribution to the four fermi interaction arising from the 
massless Higgs particle exchange. Therefore in the limit a; ^ — oo we expect to observe an 
exchange of a massless particle. 

The hypergeometric functions behave in the limit x — — oo 

^, , , r(c)r(6-a) , r{c)r{a-b) f, 

hm Fia,b,c,x) = ^), ' ^) f + ^) { ^ x^'' 



r(6)r(c-a) r(a)r(c-6) 



Tic) Tib - a) , Tic) Tia - b) . 



i [0) I [c — a) 1 [aj i (c — 0) 
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Hence /(a, b, x) for x — oo takes the form 



x— >— oo 2,71 

with 



1 \ (-lY'''x''+''Tab 0<a + 6<l 

hm —I^(a,b,x)-' { ' , (B12) 



_ r(i-a)r(i-&)r(a + 6) 

r(a)r(6)r(i-a-6) • ^'''''^ 

Using ()B12|) the amphtude ()4.4|) becomes in the hmit x — ^ — oo 

- (27r)t r|^,_2,, r|^,_2.. rt,3,-i-2.3 / ^a: a:-"'*-^ . (B14) 

J — oo 

Thus, we observe an exchange of a massless particle, which we identify as the Higgs-particle. 
Note that the prefactor in ()B14|) is the expected relative factor between the Yukawa couplings 
in string and field theory basis 0, |3| • 
Applying the limit for our second amplitude ()4.7j) we obtain 

3 ^ „ 
~ (27r)i nrL„-i-2., / clx x-'^''-' (B15) 



7=1 

and again we can observe a massless Higgs exchange in this limit. 
The amplitude {V\* V\ 

2 2 2 2 

The analysis for both amplitudes, ()4.4|1 and ()4.7|1 is similar, so that we will describe 
the steps for the first one and apply these later for the second amplitude. We start by 
investigating the integral K{9i, 02, 9^) and turn later to T{6i, 6*2, 6*3). 

Since in this interval the amplitude is finite even in the low energy limit, we do not 
have to subtract anything. Thus, we can send a' to zero and obtain 
1 

K= / da;— + /(-^^2,l + ^^2,x) / (1-^3,1 + ^^3,2;)]-^ . (B16) 

Let us split the integral ()B16|) by using the expression 

x(l — X) X 1 — X 
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Let us first evaluate the integral starting with the first summand of ()B16|) which is given by 
1 

K^= / Ax-[i{-e^A + v^,x) li-e^A + v^^x) i{i-e^A + v^,x)]--^ . (bis) 

Jo ^ 

Substituting e~* for x we obtain 

K,= dt [/(-^^i,l + z/i,e-*) /(-^^2,l + i^2,e-*) /(I -^3,1 + ^^3, e-*)] (B19) 





Mathematica is not able to evaluate this expression numerically since it is hard to maintain 
numerical precision for large t. Therefore we will split integral ()B19|) into the range from 
to T and from T to oo. For the computation of the first region we will use Mathematica to 
evaluate it numerically, while for the second region we replace the hypergeometric functions 
by their asymptotic behavior given in ()B1|) 

K^= f dt + 1^1, e-*) / (-^2, 1 + ^^2,6-*) 1(1-^3,1 + ^^3,6-*)]"^ 

^0 

/•oo ^ 

+ 7r3/2 j dt [{t + ln5(-^i, 1 + z/i)) (t + ln5(-^2, 1 + 1^2)) {t + In 5(1 - ^3, 1 + ^^3))]"^ ■ 
Let us assume that the two angles 9i and 62 are equal to each other 

61 = 62 = 6 — ^ 0s = — 2,6 . 

Then Ki simplifies to 

Ki= dt [I [-9, 1 + 29, e-^) I [-9, 1 + 29, e-*) I {1 + 29,-1-4:9,6-*)] ' (B20) 
Jo 

r°° 1 

+ 7r3/2 / &t{t + \n5{-9,l + 29)Y\t + \Yi5{l + 29,-1 -A9)Y^ . 
Jt 

Now we turn to the second term we get after splitting the integral. Again we substitute e~* 

for X, set, as above, 9i = 92 = 9 and obtain 

/■°° e"* -1 

K2= / dt- [P {-9, 1 + 29, e~') I {1 + 29,-1 -A9,e-')] ^ . 

Jo 1 — e 

As above we have to split this integral into two parts, where we replace the /'s by their 
asymptotic behavior 

K2= f dt (-^, 1 + e-*) / (1 + 29, -1 - 49, e"*)] (B21) 

Jo 1 — e 

+ 7r3/2 / (^t + ln5i-9,l + 29)y\t + \n5{l + 29,-1-49))'^ . 

Jt 1 - e * 

Applying the same procedure for the other sector we obtain 
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i<01<l i<02<l i<^3<l 



In this sector we obtain 

fT 



r -1 

Ki= dt (1-^,-1 + 2e,e"*) / (20 -1,3-4^, e"*)] ^ (B22) 

/oo ^ 
dt (t + ln5(l -0,-1 + 20))"^(t + ln5(20- l,3-40))"2 



and 

K2= f dt — ^— ^ [/^ (1 - 0, -1 + 20, e~*) / (20 - 1, 3 - 40, e"*)] (B23) 
Jo 1 — e 

+ 7r=^/2 / (t + ln5(l-0,-l + 20))-^(t + ln5(20-l,3-40))-2 . 

Jt 1 - e * 

The whole integral K{9) is given by the sum of Ki and K2. 

T(0l,02,03) 

Let us now analyze the massive string state contribution to T(0i,02,03), where in 
the field theory the proton decay takes place via Higgs particle mediation. Thus, in contrast 
to the numerical analysis for proton decay via a gauge boson exchange we observe a pole 
that corresponds to the Higgs exchange. In order to obtain the four-Fermi interaction term 
due to the massive string states, we need to subtract this pole before taking the low energy 
limit. 

Let us split the integral ()4.5p into two parts (again we assume that 0i = 02 = 0) 

/ da;x-"'"-^(l [/2 (-0,1 + 20, a;) 7(1 + 20,-1 -40, x)] ^ (B24) 
J —00 
^0 

+ J dx (1 - [P (-0, 1 + 20, x) / (1 + 20, -1 - 40, ; 

Now we replace a; by 1 — in the first summand and in the second by 



1 



in(i-L) (e-)«'« (1 _ e-)«'.+i ^ ' 

^ - [P (-0, 1 + 2 0, j^) 1(1 + 20,-1 -40,^)]-^ 



In(l--i) ) U e J 
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To simplify the computation, we break up both terms into two parts 

^ [P {-6,1 + 26,1-6') I {1 + 26,-1-46,1-6' 
dz- 



1 r°° 

JTl 

^^[P{-6,l + 26,jl^) l{l + 26,-l-A6,jl^)] 



z + \nS{-6, 1 + 26))'^ {z + \n6{l + 26, -1 - 46))~^ 



3 /"°° 

+ 772/ dz 



Here we replaced the hypergeometric expressions by their respective limits in the range 
from Ti to 00 and T2 to 00 . As mentioned above in order to get the four- Fermi interaction 
contribution, we need to subtract the ^ pole and take the low energy limit 

T(6) - linn i T H. 1 + 2^, 1 " e') I (1 + 26, -1 -46,1- 6^ 

^ r^^ [I' {-e, 1 + 26, j^)i{i + 26, -1 - A6, ^)] 



3 

+ 7r2 



In(l-i) (e-)a'«(l_e^)a't 

^ + In 1 + 2^))"^ + In 5{1 + 26, -1 - 4^))"^ 



L ' 

COO 



{z + In 5{-6, 1 + 26))-' {z + In 5{1 + 26, -1 - 4g))~2 
4 (e^)a'«(l_e-)"'* / ■ 

For the second region ^ < 61 < 1, ^ < 62 < 1 and | < 6*3 < 1, T{6) takes the form 
T{6) = hm I r - ^' -1 + ^ - ^ + ^ - ^ - 



+ (^(27r)i r,_,,_i+2. r!,^,,,3_,, d;^ (e^j-'^^+i (1 - e^)-'-^ - 
^ [I' (1 - ^, -1 + 26, ^) / (-1 + 2^, 3 - 46, ^)] 



In(l-i) (e^)a'« (1 _ e-)a't 

^ + ln5(l -6,-1 + 26))-^ {z + In 6{-l + 26,3- 46))'^ 



00 



+ n2 I dz / ^ , . s I ■ 

Jt2 (e^)" " (1 - e^)" * J 

Mathematica is not able to take that limit, however by plugging in different small values 
for a' (keep in mind that the Mandelstam variables s, t and u have to satisfy momentum 
conservation s + t + u — 0) we get a stable contribution for T{6). 
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The amplitude {V}^ V}lV}^i V}l) 



The analysis is simpler for {V_i V_i V i V i) because of the symmetry of the am- 

2 2 2 2 

plitude: after splitting the integral ()B17|) both parts give the same contribution, so that we 
only need to focus on one part and multiply by a factor of two. Following the same steps 
as above the integral M becomes 

M = 2 dt JJ Vsin[7r(l + z/7)] L'^l + Ui) + tt^^^^ / dt JJ (t + ln5(l + z//, 1 + z//))-^ 
-^0 1=1 

Replacing uj by dj and assuming that 6i = 62 we get for 

• -i<0i<O -|<e2<0 |<e3<l 

M = 2 [ dt sin[7r(l + 29)]^/sm[n{-l - 49)] L-\l + 29) L'^-l - 49) (B26) 
Jo 

+ 7r^/2 j &t{t + In (5(1 + 29,1 + 29))'^ {t + In 5{-l -49,-1- 49))~^ , 

and for 

• |<^i<l 1<92<1 l<9:i<l 

rT 

M = 2 dt sm[7T{29 - l)]^sm[n{3 - 49)] L-\29 - 1) L"5(3 _ 4^) (B27) 
Jo 

POD 

+ 7r3/2 / dt {t + \n6{29 - 1,29 - l))'\t + \n5{3 - 49,3 -49))'^ . 
Jt 
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